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USE  OF  FLOYD'S  ALGORITHM  TO  FIND  SHORTEST  RESTRICTED  PATHS 
VICTOR  KLEE  and  DAVID  LARMAN 

University  of  Washington,  Seattle  and  University  College,  London 

Abstract.  In  a directed  network  with  no  negative  circuit,  Floyd's  algorithm 
finds,  for  each  pair  of  nodes  x and  y,  a shortest  path  from  x to  y. 

Here  the  procedure  is  extended  to  minimize  more  general  length-functions  over  sets 
of  paths  that  are  restricted  in  various  ways. 

Introduction. 

Throughout  this  paper,  G denotes  a complete  directed  graph  with  n nodes  and 

2 

n - n edges.  Each  edge  is  an  ordered  pair  (i,  j)  of  nodes  and  has  as  its  length 
a number  x(i,  j)eR  = ]-«>,  «].  For  notational  convenience  G's  node-set  is 
assumed  to  be  the  set  N = {1,  ••*,  n}. 

A walk  from  x ^ y is  a node-sequence  (x^,***,  x^)  such  that  x^  = x, 

X.  = y,  and  t > 0.  It  is  a chain  if  no  node  is  repeated  and  a circuit  if  x^.  = x^ 

X — — — 1,  Q 

but  there  is  otherwise  no  repetition.  Both  chains  and  circuits  are  called  paths , a 
practice  that  is  unusual  but  is  convenient  for  our  purposes. 

Floyd's  algorithm  [R]LW][F][H2][L]  initializes  S[i,  j]  ■<-  X(i,  j) 
for  all  i,  jcN  and  then  proceeds  as  follows: 

for  k 1 until  n ^ 
for  i -f-  1 until  n do 
for  j 1 until  n ^ 

S[i,  j]  min{S[i,  j],  S[i , k]  + S[k,  j]}. 


If  there  are  no  circuits  of  negative  length  then  S[x,  y]  anerges  as  the  length  of 
a shortest  path  from  x to  y.  The  computation  is  easily  modified  to  find  shortest 
paths  in  addition  to  their  lengths. 


In  the  present  paper  the  procedure  is  extended  to  deal  with  a family  F of 
sets  of  walks  and  with  a walk-length  function  L more  general  than  the  usual  one. 
Under  suitable  assumptions  the  extended  procedure  finds,  for  each  choice  of 
X,  yeN  and  ZeF,  a shortest  Z-path  from  x to  y.  That  is,  L is  minimized 
over  the  set  of  all  paths  from  x to  y that  belong  to  Z.  In  the  "classical" 
case,  Z is  the  set  of  all  paths  (or  walks)  in  G,  F = {Z},  and  the  length  of  a 
walk  is  the  sum  of  the  length  of  its  edges. 


2.1 


The  Assumptions. 

The  function  L is  used  to  measure  the  length  of  a walk  in  terms  of  the  lengths 

★ 

of  its  edges.  It  is  assumed  the  range  of  L is  contained  in  R , the  domain  of 

★ 

L is  the  set  of  all  finite  sequences  in  R , and  the  following  two  conditions 
are  satisfied: 

■k 

(1)  ^ a-i,**',  ot^eR  and  0<s<t  then 

01^)  = L(L(q(^,-",  a^),  cx^)) 

(2)  if  ’ ^2  ^ *■  “l  ~ ^1  “2  “ ^2  1(01^,02)  ^ L(e^,62)- 

The  length  of  a walk  W = (x  ,•  • • , x ) is  defined  as 

0 W 

Lj^(W)  = L(x(x^»  x^),***,  x^))* 

By  (1),  ~ ^ ^ X to  y and  V is  a 

walk  from  y to  z.  Here  UV  denotes  the  walk  that  follows  U from  x to  y and 

then  follows  V from  y to  2. 

Among  the  admissible  functions  L are 

Lp(a^,-*-,  a^)  = (a^+* • for  an  integer  p > 0, 

*-oc.(“l  *' ’ * • “t^  ~ max(a^,'**,  a^), 

L^(a^,*’*,  a^)  = da-i  |^+’*'  + |oi^|'’)^'^^  for  a real  p > 0. 

The  usual  L is  L-j.  The  function  is  also  of  practical  interest,  for  if  G 

is  initially  equipped  with  nonnegative  real  edge-weights  Y(i»  j)  representing 
flow  capacities  and  if  x(i,  j)  = -Y(i.  j)  for  all  i and  j,  then  the  shortest 
paths  with  respect  to  L are  those  of  maximum  flow  capacity  for  specified  initial 

00 

and  terminal  nodes  [H1][H2][L]. 


2.2 


When  y and  V are  sets  of  walks  let 

yy  = {UV:  Uey,  VeV,  U ends  where  V starts}. 

Thus  a walk  x^)  belongs  to  yv  if  and  only  if  there  exists  s such  tliat 

0<s<t,  (Xq,-'‘,  x^)eU  and  (x^,---,  x^)cy.  The  first  assumption  about  F is: 

(3)  vf  0<s<t  and  (x^,**-,  x^)eWcF  then  there  exist  y,  V such  that 
(Xo»--*,  Xg)eyeE.  (Xg,-*-,  x^)ey£F,  and  yV  c W. 

For  a simple  but  interesting  example,  suppose  that  a set  of  special  edges 

of  G is  given,  and  for  0<k<£  let  W(k)  denote  the  set  of  all  walks  (Xq,--*,x^) 

such  that  is  special  for  at  most  k values  of  i.  Let  E = 

{W(k)  = 0<k<4}.  Note  that  W(i)W(j)  c W(k)  when  i+j<k. 

The  edges  of  a walk  W = (x^,- • • , x^)  are  (x^,  x^ ),•••,  (x^_1 , x^) . A path 

P is  associated  with  W if  P also  starts  at  x and  ends  at  x^,  P is  a 
0 t 

subsequence  of  W,  and  each  edge  of  P is  an  edge  of  W.  It  is  assumed  E»  L 
and  A are  interrelated  as  follows: 

(4)  if  WeWcF  and  P is  a path  associated  with  W then  PcW  and  L,(P)<L,(W). 

Since  there  are  only  finitely  many  paths  in  G,  a consequence  of  (4)  is: 

(5)  For  each  x,  yeN  and  ZeE,  either  there  is  no  Z-walk  from  x to  y or 
there  is  a Z-path  which  is  a shortest  Z-walk  from  x to  y. 

As  is  shown  in  the  last  section  of  the  paper,  a number  of  problems  on  shortest 
restricted  paths  can  be  formulated  and  efficiently  solved  in  terms  of  a function 
L and  a family  F of  sets  of  walks  satisfying  conditions  (1)  - (4). 


r 


^ 

3.1 

The  extended  version  of  Floyd's  algorithm  (EVFA)  starts  with  the  n^n  matrix 
A of  edge-lengths  of  G,  procedures  for  computing  L(a)  and  L(a,  S)  for  all 
ct,  BeR  , and  a suitable  representation  of  the  family  F of  sets  of  walks.  Also 
required  is  a set  T of  triples  (U,  V,  W)  of  members  of  F such  that: 

(6)  UV c W for  each  (U,  V,  W)cT; 

I 

I 

(7)  if  0<s<t  and  (x  x.)€WeF  then  there  exist  U and  V such  that  j 

(x  ,•••,  X JeU,  (x  ,•••,  X JeV,  ^ (U,  y,  W)£T. 

U 3 *"  i L *“  “ 

Of  course  T may  be  taken  as  the  set  of  all  triples  (U,  V,  W)  of  members  of  F | 

such  that  yy  c W,  but  it  is  most  efficient  to  have  T as  small  as  possible  subject 

to  (6)  and  (7).  (In  the  example  following  (3)  in  the  preceding  section,  it  would  j 

be  best  to  let  I = {(W(i),  W(j),  W(k)):  i+j  = k<£  } rather  than  using  i+jsk.) 

★ 

For  each  ZeF  the  algorithm  outputs  four  nxn  matrices:  an  R -valued  S^, 
an  integer-valued  an  F-valued  and  an  F-valued  V^-  At  the  time  of  out- 

put these  satisfy  the  following  conditions  for  all  x,  yeN: 

(8)  if  there  is  no  Z-path  from  x y then  S^Cx,  y]  = ” and  M^Cx,  y]  = -1; 

(9)  if  there  is  a Z-path  from  x y then 

(a)  S^Lx,  y]  (possibly  «)  length  of  a shortest  Z-patli  from  x ^ y, 

(b)  M^Cx,  y]  = 0 vf  (x,  y)  is  such  a shortest  path;  otherwise,  M^Cx,  y]  is 
the  i ndex  m of  ^ intermediate  node  on  such  £ shortest  path , and  the  path 
itself  is  formed  from  a shortest  U,[x,  y]-path  from  x to  m followed  by 
a shortest  V^Cx,  y]-path  from  m to  y. 

Using  the  output  M,  U and  V of  EVFA  the  path-tracing  algorithm  (PTA)  actually 
finds  the  shortest  paths. 


3.2 


EVFA:  EXTENDED  VERSION  OF  FLOYD'S  ALGORITHM 


begin 

for  i ■<-  1 until  n do 
for  j 1 until  n do 
for  each  WeF  ^ 
if  (i.  j)eW  then  begin 


S^[i,  j]  L(x(i,  j)); 

M^[i,  j]  0 


else  begin 


Sy[i , j]  ^ 

M^j[i,  j]  ^ -1 
end  of  initialization; 


for  k -f-  1 until  n do 


for  i 1 until  n do 
for  j -t-  1 until  n do 
for  each  (U,  V,  W)€T  ^ 

^ (L(S^j[i,  k],  Sy[k,  jj)  < Sy[i,  j]) 


or  (M|^[i,  j]  = -1  ^ M^j[i,  k]  ^ -1  and  My[k,  j]  ^ -1) 
then  begin 


Sj^[i,  .i]  L(S^[i,  k],  S^[k,  j]); 

f\[i , j]  ^ k; 

U|j[i . j]  U; 

yi.  j]  - y 

end  of  main  loop; 


3.3 


for  each  WeF  ^ begin 


print  S^,  and 


end 


end 


t 

I 

; In  the  path-tracing  algorithm,  STACK'S  members  are  alternately  node-indices 

i and  members  of  F.  (It  is  often  convenient  in  practice  to  represent  the  members 

i of  F by  negative  integers.)  When  a shortest  Z-path  from  x to  y is  desired, 

r 

I STACK  is  initialized  as  (y,  Z,  x).  As  STACK  is  processed,  node-indices  are  added 

1 to  PATH,  which  emerges  as  a shortest  Z-path  from  x to  y.  (PTA  has  no  output 

i when  there  is  no  Z-path  from  x to  y.) 


PTA:  PATH-TRACING  ALGORITHM 


begin 

if  M^Cx,  y]  = -1  then  goto  NONE; 
STACK[1]  ^ y;  STACK[2]  ^ Z;  STACK[3]  ^ x; 
s-«-3;  p-<-0; 
while  s ^ 3 do 
begin 

i ^ STACK[s]; 

W ^ STACK[s-l]; 
j ^ STACK[s-2]; 

m ^ Mjj[i , j]; 

i_f  m = 0 then  begin 

P ^ P+1 ; 

PATH[p]  i ; 
s f s-2 


f 

) 

\ 

3.4 

i 

end 

else  begin 

STACK[s-l]  ^ Vy[1,  j]; 

STACK[s]  ^ m; 

STACK[s+l]  ^ Uy[i , j]; 
s s+2; 

STACK[s]  ^ 1 
end 

end  of  loop; 

PATH[p+l]  ^ STACKfl]; 
print  PATH; 

NONE: 

end 

In  describing  the  efficiency  of  EVFA  and  PTA  we  use  the  uniform  cost  criterion 
in  the  RAM  model  of  random  access  computation  [AHU].  Initialization  of  EVFA  requires 
time  0(|F|n  ) and  each  passage  through  the  main  loop  requires  time  0(|T|n  ), 

SO  the  overall  time-complexity  of  EVFA  is 

0(lF|n^  + iTln^). 


This  counts  each  evaluation  of  L(a)  or  L(oi,  e)  as  a single  step. 

For  each  choice  of  x,  y,  Z,  PTA  requires  time  0(n)  to  find  a shortest 
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Z-path  from  x to  y.  Thus,  starting  from  the  output  of  EVFA,  time  0(|F|n  ) 
is  required  to  find,  for  all  x,  yeN  and  ZcF,  a shortest  Z-path  from  x to  y. 


It  remains  to  show  that  EVFA  and  PTA  compute  what  is  claimed  for  them.  In 
the  case  of  PTA,  this  follows  from  a routine  inductive  argument  showing  that  at 


3.5 


the  end  of  the  initialization  and  also  at  the  end  of  each  passage  through  the 
loop  there  exists  a shortest  Z-path  P from  x to  y such  that  the  following 
three  conditions  are  satisfied: 


(a)  s is  odd;  STACK[1]  = y;  if  p > 0 then  PATH[1]  = x;  i 

(b)  PATH[1 ],•••,  PATH[p],  STACK[s],  STACK[s-2],---,  STACK[1]  is  a subsequence  ! 

I 

of  P;  I 

i 

(c)  for  each  odd  k with  3<k<s,  the  segment  of  P that  joins  i = STACK[k]  ' 

j to  j = STACK[k-2]  is  a shortest  STACK[k-l ]-path  from  i to  j. 

We  turn  now  to  EVFA.  For  WeF  and  0<k<n,  let  Wj^  denote  the  set  of  all 

walks  (x  x.)eW  such  that  x <k  when  0<s<t.  Thus  the  end  nodes  of  walks 

in  W|^  are  unrestricted  but  all  intermediate  nodes  are  in  {I,***,  k}.  Let 

F.  = {W.  : WfiF}.  Since  W = W for  all  WeF,  it  suffices  to  prove  the  following 
sK  - K n = 

for  0<k<n: 

(10|^)  After  the  k*"  passage  through  the  main  loop  of  EVFA,  condi tions  (8)  and 
(9)  are  satisfied  for  all  x,  yeN  and  ZeFj^. 

The  proof  is  by  induction  on  k.  Here  initialization  is  regarded  as  the  0^*^  passage 
through  the  main  loop  and  assertion  (10^^)  is  obvious  because  is  merely  the 
set  of  all  edges  (paths  (x^,  x^))  in  Z. 

From  the  argument  below  it  follows  that  {10|^)  holds  for  all  k regardless  of 
the  order  in  which  i,  j and  (U,  V,  W)  appear  in  the  main  loop.  In  particular, 
the  main  loop  could  be  written  as 

i J 

for  k t-  1 until  n ^ j 

for  each  (U,  V,  W)eT  ^ ] 

for  each  (i,  j)e{l,**',  n}x{l,'*-,  n}  ^ . . ..  I 


That  is  convenient  for  programming  in  some  languages  with  special  array-handling 
capabilities,  such  as  APL. 


3.6 

Now  suppose,  with  0<k<n,  that  holds,  and  consider  the  k ^ passage 

through  the  main  loop.  We  note  first  that  if  i or  j is  k then  there  do  not 
exist  U,  VeF  such  that  UV  c W and 

(d)  M|^[i,  j]  = -1  but  M^[i,  k]  -1  My[k,  j]  or 

(e)  L(Sjj[i,  k],  Sy[k,  j])  < S^^Ci,  j]. 

Suppose,  for  example,  that  j = k.  If  (d)  holds  there  is  no  W|^_i-path  from  i to 
k but  there  is  a U|^_i-path  U from  i to  k and  there  is  a y|^_i-P3th  V from 
k to  k.  But  then  UV  is  a Wj^-walk  from  i to  k and  by  condition  (4)  there 
is  an  associated  Wj^-path  P from  i to  k.  Plainly  P is,  in  fact,  a W|^_.j-path, 
and  that  is  a contradiction.  A similar  contradiction  is  derived  from  (e),  using 
conditions  (1)  - (4)  and  the  inductive  hypothesis.  It  follows  that  the  k ^ rows 
and  k*^^  columns  of  and  are  unchanged  by  the  k^^  passage  through 

the  main  loop.  Hence  (iOj^)  holds  for  all  and  x,  yeN  with  x = k or  y = k. 

The  case  in  which  remains. 

Supposing,  still,  that  0<k<n  and  (lOj^  .| ) holds,  consider  WeF  and  i,  jeN 
with  ij^k^^j.  We  discuss  only  the  case  in  which  (e)  holds  at  some  time  during  the 
k ^ passage,  for  the  other  cases  (described  in  terms  of  (d)  and  (e))  are  similar.  Let 

u = min{L(Sjj[i,  k],  Sy[k,  ,1]):  (U,  V,  W)eT}, 

^ u 

the  minimum  during  the  k passage,  and  let  (U',  V')  be  the  first  pair  (U,  V) 
for  which  the  minimum  is  attained.  Then  at  the  end  of  the  k^*^  passage, 

Sy[i,j]  = UiM^^[i,  j]  = k,  Uy[i,  j]  = y,  i,  j]  = V'  . 

Let  U be  a shortest  u!,  i-path  from  i to  k and  let  V be  a shortest  V, 
path  from  k to  j,  whence  L (UV)  = y.  Then  UV  e U'V  c W and  hence  UV  is 


I a W|^-walk  from  i to  i.  For  (lOj^)  it  suffices  to  show  UV  is  a shortest 

W|^-path  from  i to  j.  Consider  an  arbitrary  shortest  Wj^-walk  W from  i to 
1 j and  an  arbitrary  associated  path  P = (Xq,---,  x^).  Then  P e and 

L.(P)  ^ L,(W),  whence  of  course 

A A 

' 4(P)  = L^(W)  ^ L^(UV). 

j 

The  node  k appears  in  P for  otherwise  it  is  true  at  the  end  of  the  (k-1)'^" 

I th 

I passage  that  S^Ci>  j]  = never  holds  during  the  k " passage, 

contrary  to  hypothesis.  With  = k there  exist  U and  V such  that 
[ (x  X ) e U,  (x  X ) e y and  (U,  V,  W)  e T.  But  then 

I 

I L^(UV)  < L(Sy[i,  k],  Sy[k,  j]  < L(L^(x^.-.-,  x^),  ^(x^,---,  x.J)  = L(P), 

; whence  L (UV)  = L (W)  and  UV  is  a shortest  W -walk  from  i to  j.  If  UV  is 

1 not  a path  it  has  an  associated  path  that  misses  k,  and  that  was  shown  to  be 

i 

j impossible. 

i 

j 

i 

r 

[ 

I 

I 

j 

I 

i 

i 

j 

i 

t 

! 

! 
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Though  conditions  (l)-(4)  suffice  for  the  validity  of  EVFA,  some  additional 
conditions  aid  in  verifying  condition  (4)  for  specific  applications.  The  function 
L is  said  to  be  nice  if  in  addition  to  (1)  and  (2)  it  satisfies  the  following  two 
conditions: 

(11)  each  point  of  rng  L is  fixed  under  L;  that  is,  LCa-j."',  a^)  = 

L(L(oii » “ ■ , a^) ) ; 

(12)  rf  3^0  then  L(a,  6)  ^ L(oi)  < L(3>  a)  f^  £0^  a e rng  L. 

Note  that  each  of  L„,  L and  is  nice. 

P oo 

If  W is  a walk  (Xq.'”.  x^)  and  the  proper  segment  (x^,---,  x^)  of  W 
is  a circuit  C then  C is  called  an  intermediate  circuit  of  W and  denotes 

the  walk  that  remains  when  all  of  C but  x^  or  is  removed  from  W.  More 
precisely,  when  0 < r is  the  walk  (x^,-*-,  x^,  x^),  (x^,---,  x^,x^^^) 

or  (x^,---,  x^)  according  as  s+l<t,  s+1  = t or  s = t,  and  when  s < t 

is  the  walk  (x^,---,  x^_^ , x^,---,  x^) , (x^_i,  x^,---,  x^)  or  (x^,---,  x^)  accord- 

ing as  r-l>0,  r-1  = 0 or  r = 0. 

Note  that  a walk  is  a path  if  and  only  if  it  has  no  intermediate  circuit.  Hence 
condition  (4)  can  be  deduced  from  repeated  application  of  the  following  condition: 

(13)  W = {Xq»---,  x^)  e W e F and  (x^,--*,  x^)  is  an  intermediate  circuit 

of  W then  W e W and  L,(W  ) < L, (W). 

— rs  ~ X'  rs  X 

Note  that  the  inequality  of  (13)  always  holds  when  L = L^.  In  other  cases  it 

can  often  be  deduced  from  the  following  result. 

(14)  If.  L is  nice,  C=(x^,-",x)  iim  intermediate  circuit  of  a^  walk 

W = (Xp,'”,  x^)  ^ L^(C)  s 0 then  L;^(Wy.5)  ^ L^(W). 

To  prove  (14),  note  that  W = CW  if  r = 0,  W = W C if  s = t,  and  if 

I 5 I 

0<r<s<t  there  are  walks  U and  V such  that  W = UCV  and  W . = UV.  We  consider 

rs 


only  the  third  case  for  the  others  are  similar  to  it.  With  L^(C)  > 0,  it  follows 
from  (11),  (12)  and  (1)  that 


L^(U)  = L(L^(U))  < L(Lj^(U).  L^(C))  = L^(UC), 


then  from  (1)  and  (2)  that 

4(W^s^  =4(UV)  = L(L^(U),  L^(V))  ^ L(L^(UC),  L^(V))  = L^(UCV)  = L^(W). 

Below  are  some  illustrative  problems  on  shortest  restricted  paths.  In  each 
case,  "find  shortest  paths"  means  that  for  each  x,  yeN,  either  a shortest  path 
from  X to  y (among  those  satisfying  the  indicated  restrictions)  must  be  found 
or  it  must  be  concluded  that  no  path  from  x to  y satisfies  the  restrictions. 

(A)  A set  of  edges  is  given.  For  each  k < £,  find  shortest  paths  that  use 

at  most  k of  the  special  edges. 

(B)  A set  of  nodes  is  given.  For  each  k < £,  find  shortest  paths  that  use  at 

most  k of  the  special  nodes. 


(C)  A sequence  of  s sets  is  given,  each  consisting  of  nodes  or  edges  or  a 


mixture.  For  each  choice  of  (k,  ,•••,  k ) with  k < I for  all  r,  find 

I s r r 


shortest  paths  that  use  (for  all  r)  at  most  k^  of  the  elements  of  the  r^*^  set. 


(D)  In  addition  to  the  R -valued  edge-lengths  x(x,  y),  integer  edge-lengths  TT(x,y) 


^ 0 are  given.  Each  walk  has  its  usual  length  L and  also  a length  I 

A "I 


where  I is  L . An  integer  i > 0 is  given.  For  each  k < £,  find  L -shortest 
1 ^ 


paths  P subject  to  the  restriction  that  I.(P)  s k. 

A 


(E)  The  nodes  of  G are  partitioned  into  two  disjoint  sets  A and  B,  and  an 


integer  1^0  is  given.  For  each  k $ find  shortest  paths  that  oscillate  at 
most  k times  between  A and  B. 
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(F)  A subgraph  H of  G and  an  integer  I > 0 are  given.  For  each  k < £, 
find  shortest  paths  P for  which  P n H has  at  most  k components. 

(G)  A set  M of  edges  of  G is  given.  Find  shortest  M-alternating  paths. 

As  can  be  seen  by  reference  to  conditions  (1)-(4)  and  (13)-(14),  the  discussions 
of  (A)-(G)  below  are  valid  (that  is , EVFA  can  be  applied  for  the  stated  purpose) 
rf  L i^  and  also  if  L is  nice  and  L^(C)  > 0 for  each  circuit  C inter- 
mediate to  ^ walk  belonging  to  ^ member  of  F , where  F the  family  of  sets  of 
walks  used  for  the  particular  problem.  (Problem  (G)  requires  an  additional  condition, 
stated  later. ) 

(A)  This  problem,  which  was  mentioned  earlier,  is  straightforward.  For  0<k<Z, 
let  W(k)  denote  the  set  of  all  walks  (x^,--*,  x^)  such  that  the  edge 
is  special  for  at  most  k values  of  i.  Let  F = {W(k):  0<k<£}.  Let  T = 
where 


T,^  = {(W(i).  W(j).  W(k)):  i+j  = k}. 

Then  |I|^|  = k+1  and  jT]  = (£+1 ) (£+2)/2.  The  overall  time-complexity  of  EFWA  for 

2 

this  problem  is  0(£  n ). 

(B)  This  problem  is  similar  to  (A),  but  it  is  included  to  illustrate  the  way 
in  which  the  end  behavior  of  walks  must  sometimes  be  considered  in  constructing 
F and  T for  the  application  of  EFWA.  For  0<k<£  let  W__(k)  (resp.  W_^(k), 
W^_(k).  W^^.(k))  denote  the  set  of  all  walks  such  that  the  node  x. 

isspecial  for  at  most  k values  of  i with  0<i<t  and,  in  addition,  neither  x^ 
x^  is  special  (resp.  x^  is  special  but  x^  is  not,  x^  is  special  but  x^ 
is  not,  both  x^  and  x^  are  special).  Let  F consist  of  the  sets  W__(k)  for 
k < £ , the  sets  W_^(k)  and  W^_(k)  for  k<£-l , and  the  sets  W^^(k)  for 
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i 


t 

I 


k < t - 2.  Let  T = l/l,  , where  I,  consists  of  the  triples 
= o=k  -K 

(W— (i)>  W__(k))  for  i+j  = k, 

W_^(j).  W_+(I<-1)).  (W+.(i),  W_.(j).  W^_(k-1)) 

and  (W_+(i),  W^.(j),  W__(k))  for  i+j  = k-1 , 

(W_+(i).  W++(j).  W_+(k-l)),  (W^+(i).  W+.(j).  W^_(k-1)) 

and  (W^_('i).  W^^(k-2))  for  i+j  = k-2, 

and  (W++(1).  W++(j)»  W++(k-2))  for  i+j  = k-3. 

Again,  |T|  is  0{t^)  and  the  complexity  of  EVFA  is  0(£^n^). 

(C)  This  is  included  to  illustrate  the  application  of  EVFA  when  the  desired 

a. 

paths  are  subject  to  several  restrictions.  In  order  to  avoid^notational  morass, 
only  the  case  of  sets  of  edges  is  discussed.  For  k^  < k^  < let 

W(k^,*--,  k^)  denote  the  set  of  all  walks  (x^,-*-,  x^)  such  that,  for  l<r<s, 
(Xi_i,  x^. ) belongs  to  the  r set  of  edges  for  at  most  k^  values  of  i.  Let 
T consist  of  all  triples 

ij).  W(j^,---,  j^),  W(k^,---,  k^)) 

such  that  for  lsr<s,  i + j = k < t . Then 

r ^r  r r 

III  ' Vo  Jk2=0  ••• 

and  the  complexity  of  EVFA  is 

(D)  This  may  be  regarded  as  the  integer-weighted  version  of  a problem  of  which 
(A)  is  the  cardinality-weighted  version.  Similar  extensions  are  available  for  the 
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other  problems  considered  here.  For  Osk<£,  «let  W(k)  be  the  set  of  all  walks 

W for  which  I (W)  < k.  Define  F and  T in  the  obvious  ways.  The  complexity 
2 2 

of  EVFA  is  0(^-  n ).  For  a closely  related  treatment  of  this  problem  and  of 
(A),  see  the  discussion  of  the  Bellman-Ford  method  in  [L,  pp.  74-75,  92-931. 

EVFA  is  similar  to  the  Bellman-Ford  method  but  is  more  general.  Roughly  speaking, 
it  amounts  to  replacing  the  additive  semigroup  {0,1,2,...}  of  Bellman-Ford 
by  an  arbitrary  semigroup. 

(E)  This  is  a special  case  of  a more  general  problem,  which  may  be  fonnulated 
as  follows:  A function  <))  is  defined  on  a set  of  nodes  and  edges  of  G,  with 
rng  c {!,•••,  ml,  and  an  integer  I > 0 is  given.  For  each  k < £,  find 
shortest  paths  along  which  ((>  has  at  most  k relative  extrema. 

As  the  term  is  used  here,  a relative  extremum  of  a real  sequence  a^) 

is  an  ordered  pair  (r,  s)  such  that  0<rss<u  and 

r-1  r s s+1  r-1  r s s+1 

For  a walk  W = (x^,---,  x^),  let 

exp  W = ^^o’  ^1^’  ^1 ’ " ' " ’ ^^t-1 ’^t^ ’ ^t^ ’ 

the  expanded  version  of  W in  which  nodes  and  edges  alternate.  Let  denote 
the  sequence  of  4>-values  corresponding  to  the  elements  of  exp  W that  belong  to 

dmn  ii,  and  let  p (W)  denote  the  number  of  relative  extrema  of  W . The  general 

<J>  <}> 

problem  is  to  find  shortest  paths  P for  which  p^(P)  ^ k.  Problem  (E)  is  the 

<5> 

special  case  in  which (*) 

(*)  m = 2 and  dmn  (}>  = N = A u B,  with  <()  = 1 on  A and  1 = 2 on  B. 
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As  is  shown  below,  this  can  be  handled  by  EVFA.  However,  we  do  not  know  how  to  use 
EVFA  efficiently  for  the  general  problem,  or  even  for  the  following  special  cases: 

m = 3 and  dmn  = N; 

m = 2 and  dmn  is  a proper  subset  of  N; 
m = 2 and  dmn  (f  is  the  set  E of  all  edges  of  G. 

In  each  case  there  is  difficulty,  even  when  L is  and  all  values  of  the  edge- 
length  X are  positive,  in  constructing  a suitable  family  F satisfying  conditions 
(3)  and  (4). 

Now  let  us  return  to  (E)  in  the  formulation  provided  by  (*),  except  that  the 

condition  dmn  (}>  = N may  be  replaced  by  dmn  ()i  = N.  For  0<k<£  and  u,  ve{l,2}, 

let  W (k)  denote  the  set  of  all  walks  W = (x„,-*-,  x.)  such  that  (j)(x  ) = 

-•UV  0 t 0 

u,  <)>(k.)  = V,  and  W has  at  most  k relative  extrema  (equivalently,  W oscillates 
t (j) 

at  most  k times  between  A and  B).  Let 

F = {W„  ^(k):  u,  V£{1,  2},  k < 1} 
and  let  T consist  of  all  triples. 

^-uu^^^’-uu^'^^’-uu^*^^^  ue{l,2}  and  ({i,j}  = {0,k}  or  (i>0<j  and  i+j  = k-1)), 

(yuv^’^’-vv^'^^’-uv^'^^^  ^ and((i,j)  = (k,0)  or  (j>0  and  i+j  = k-l)X 

^-uu^^^’-uv^’^^’-uv^'^^^  ^ and((i,j)  = (0,k)  or  (i>0  and  i+j  = k-1)). 

2 3 

Then  EVFA  can  be  applied,  solving  problem  (E)  in  time  0(^  n ). 

(F)  Define  <()  on  all  nodes  and  edges  of  G,  with  <f>  = ] on  nodes  and  edges  of 
the  graph  H and  (j)  = 2 otherwise.  With  the  as  in  the  preceding  paragraph, 

the  paths  P for  which  P n H has  at  most  k components  are  precisely  the  paths  in 

L, _____ 
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W^^(2k-3)  u W^2(2k-2)  u W2i(2k-2)  u W22(2k-3). 

2 3 

Hence  (F)  can  also  be  handled  by  EVFA  in  time  0(£  n ). 

(G)  This  problem  is  also  discussed  in  a more  general  setting.  With  <}>  as  in  the 
discussion  of  (E)  and  with  m = 2,  let  a walk  W be  called  ^-alternating  if  the 
sequence  W alternates  between  1 and  2.  How  can  shortest  ^-alternating  paths  be 

found?  Problem  (G)  is  the  special  case  in  which  dmn  ij)  = E,  = 1 on  M,  and 
(P  =2  on  E ' M. 

When  dmn  <{)  =5  N,  EVFA  can  be  applied  by  taking  F = {W.|^  ,^^2*^21  >^22^  ’ where 
W^^  is  the  set  of  all  cji-alternating  walks  W such  that  the  sequence  W^  starts 
with  u and  ends  with  v.  Then  let  T consist  of  all  triples  (^uv’^vu’-uu^ 
^-uv’-vv’-uv^  for  all  u,ve{l,2}. 

Now  consider  the  case  in  which  dmn  ij)  = E.  Then  each  pair  (i,  j)  forms  a 
(j)- alternating  path  (recall  the  standing  hypothesis  that  G is  the  complete  graph 
on  N)  but  of  course  we  are  interested  only  in  paths  of  finite  length.  Define  F 
by  restricting  the  of  the  preceding  paragraph  to  include  only  walks  of 

finite  length,  and  assume 

(t)  each  alternating  circuit  of  finite  length  has  an  even  number  of  edges. 

Then  problem  (F)  can  be  handled  by  EVFA  with  T consisting  of  all  triples 

(Wup»  Wqy*  i'or  u,  v e {0,1}  and  {p,q}  = {0,1}.  However,  this  approach 

may  fail  when  (t)  fails  for  then  a path  associated  with  an  alternating  walk  need 
an 

not  be^alternating  path  and  thus  condition  (4)  may  fail.  For  example,  consider 
Fig.  1 and  note  that  (x2.X2»X^,Xg)  is  an  alternating  circuit  according  to  our 
definition,  where  the  solid  edges  are  those  in  M. 


I 
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Fig.  1:  The  walk  x^)  is  alternating  but  the  associated  path 

(Xq,x^ ,x^)  is  not. 

Brown  [B]  suggests  a method  for  finding  shortest  M-alternating  paths  in  a directed 
graph  D = (N,  E)  (no  longer  assumed  complete).  Another  directed  graph  D*  is 
constructed,  having  two  nodes  x'  and  x"  for  each  node  x of  D,  and  the  edges 
of  D are  obtained  as  follows  for  each  edge  (x,  y)  of  D with  length  x(x,y)<“>. 

•ff 

when  (x,y)  e M,  (x',y")  is  an  edge  of  D with  length  a(x,  y); 
when  (x,y)  it  M,  (x",  y')  is  an  edge  of  D with  length  a(x,  y). 

Itis  claimed  [B]  there  is  a natural  one-to-one  correspondence  between  alternating  paths 

* 

in  D and  ordinary  paths  in  D . Thus  the  problem  of  finding  shortest  alternating 

paths  in  D is  equivalent  to  the  problem  of  finding  shortest  ordinary  paths  in 

D*.  The  claim  is  correct  when  (+)  holds  but  not  in  general,  as  can  be  seen  from 

Fig.  1.  For  example,  if  D is  the  graph  of  Fig.  1 then  the  path  (x|^ ,x|  ,X2,X2 ,X2,x!j' ,x^ ' 
★ 

in  D corresponds  to  the  walk  (x^,Xi  in  D. 

In  general.  Brown's  construction  does  produce  a one-to-one  correspondence  between 
★ 

the  walks  in  D and  the  alternating  walks  in  D.  If  D has  no  negative  alternating 
★ 

circuit  then  D has  no  negative  circuit  and  the  Floyd-Warshal 1 algorithm  can  be 

★ 

applied  to  find  shortest  paths  (=shortest  walks)  in  D and  hence  alternating  walks 
in  D.  The  latter  may  or  may  not  be  paths.  For  general  graphs,  even  when  L = L-] 
and  may  or  may  not  be  paths.  For  general  graphs,  even  when  L = Li  and  X > 0,  we 


ft 

jt- 


do  not  know  how  to  apply  EFVA  directly  to  find  shortest  alternating  paths. 
However,  the  more  complicated  "blossom"  methods  of  Edmonds  [E1][E2]LL]  will 
apparently  apply  to  this  problem. 

We  close  with  a query.  For  0<!l<m<n  and  for  p e {<,  =,  >}  let 
P^(£,m,n,p)  {resp.  P^(£,m,n,p))  denote  the  following  problem: 

A complete  graph  G is  given,  with  n nodes  and  positive  edge-lengths. 

In  addition,  a set  of  m special  nodes  (resp.  edges)  of  G is  given.  Find 
shortest  paths  P in  G such  that  the  number  of  special  nodes  (resp.  edges) 
used  by  P is  in  the  relation  p to  £. 

We  have  seen  here  that  the  extended  version  of  Floyd's  algorithm  solves 
E^(Jl,m,n,<)  and  Ej^(Jl,m,n,<)  in  time  0(£  n ).  By  contrast,  the  problem 
P^(n,n,n,=)  is  essentially  the  traveling  salesman  problem  and  hence  is  NP- 
complete  rKlirKZirAHU].  What  else  of  interest  can  be  said  about  the  computa- 
tional complexity  of  EQ(^>ni»n,p)  and  E(Jl,m,n,p)  for  p e {=,  >}? 
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